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Limiting  of  Signals  in  Random  Noise 

PRAVIN  C.  JAIN,  MEMBER,  IEEE 


Abstract — The  effect  of  ideal  bandpass  limiting  on  a  signal  lying  in 
narrowband  Gaussian  noise  is  analyzed.  General  analytic  expressions  for 
the  limiter  output  components  are  derived  using  an  integral  representation 
for  the  limiter  characteristic.  This  method  allows  retention  of  the  phases 
of  all  the  signals  and  the  intermodulation  products  at  the  limiter  output, 
which  are  destroyed  in  the  characteristic-function  method  generally  used 
in  limiter  studies.  Expressions  for  the  desired  signals  and  intermodulation 
product  amplitudes  are  obtained  for  the  case  when  the  limiter  input 
consists  of  three  angle-modulated  sinusoids  and  noise.  The  analysis  is 
extended  to  n  modulated  sinusoids  plus  noise,  and  approximate  expressions 
for  the  signal,  intermodulation  product,  and  noise  terms  are  derived. 
Numerical  results  are  presented  for  the  signal  suppression  and  the 
limiter  output  signal  amplitudes  for  the  case  of  three  input  signals,  two 
of  equal  amplitude. 

Introduction 

HE  calculation  of  the  output  of  a  limiting  device  when 
its  input  consists  of  a  sum  of  several  signals  has  been 
the  subject  of  a  great  deal  of  theoretical  analysis  and  has 
resulted  in  a  number  of  widely  quoted  publications.  Essen¬ 
tially,  two  methods  of  approach  have  been  used  to  analyze 
the  problem.  The  first  approach,  commonly  known  as 
the  characteristic-function  method  of  Rice,  involves  com¬ 
puting  the  autocorrelation  function  of  the  limiter  output 
and  then  taking  the  Fourier  transform  to  obtain  the  power- 
density  spectrum.  Although  a  general  expression  for  the 
limiter  output  autocorrelation  function  can  be  derived,  its 
computation  becomes  extremely  involved  when  modulation 
of  the  input  signals  is  considered.  The  difficulty  lies  in 
determining  the  characteristic  function  of  the  signals  with 
arbitrary  modulation.  However,  if  the  signals  are  statisti¬ 
cally  independent  angle-modulated  sinusoids,  and  only  the 
average  power  or  the  magnitude  of  the  signal  and  cross- 
product  terms  at  the  limiter  output  is  of  interest,  the 
modulation  of  the  input  signals  can  be  ignored,  which 
considerably  simplifies  the  analysis.  Davenport  [1]  was 
first  to  use  this  approach  to  investigate  the  effect  of  hard 
limiting  a  single  sinusoidal  signal  and  narrowband  Gaussian 
noise.  Jones  [2]  used  the  same  method  to  analyze  the  case 
of  two  sinusoids  plus  noise.  More  recently,  Shaft  [3]  and 
Gyi  [4]  independently  extended  the  analyses  to  include  n 
sinusoids.  The  magnitude  of  any  signal  or  cross  product 
is  given  by  an  untabulated  infinite  integral  that  has  been 
numerically  evaluated  for  a  number  of  cases  of  interest. 

In  many  practical  applications,  the  phases  of  the  signal 
and  cross-product  terms  are  also  of  interest.  For  example, 
in  an  FM  or  PM  system,  the  spreading  of  the  power 
spectral  density  of  the  signal  and  cross-product  terms  will 
be  a  function  of  their  phase  modulation,  even  though  the 
average  power  in  any  one  output  component  is  indepen- 
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dent  of  its  phase.  Thus,  modulation  of  the  input  signals 
must  be  considered  in  the  analysis,  in  order  to  retain  the 
phase  of  the  signal  and  cross-product  terms  at  the  output 
of  the  limiter. 

The  second  method,  which  can  be  referred  to  as  the 
Fourier-expansion  approach,  is  carried  out  entirely  in  the 
time  domain,  in  contrast  to  the  autocorrelation  approach, 
which  treats  the  problem  in  the  frequency  domain.  Analysis 
in  the  time  domain  does  not  lose  the  phase  of  the  signal 
and  cross-product  terms.  It  also  allows  removal  of  the 
assumption  of  statistical  independence  between  the  input 
signals.  This  approach  was  used  by  Granlund  [5]  and  later 
by  Baghdady  [6]  to  investigate  the  noiseless  case  of  two 
sinusoids  passed  through  an  ideal  bandpass  limiter.  More 
recently,  Sollfrey  [7]  used  the  same  approach  to  analyze 
the  effect  of  hard  limiting  on  a  sum  of  three  or  four  sinus¬ 
oidal  signals  without  noise.  Closed-form  analytic  expressions 
for  the  amplitude  of  the  desired  signal  terms  were  obtained 
for  three  input  signals,  two  of  equal  amplitude,  and  for 
four  signals,  having  two  pairs  of  equal  amplitude.  The 
Fourier  expansion  method,  however,  has  the  drawback 
that  it  is  difficult  to  consider  the  effect  of  noise  present 
at  the  input  to  the  limiter.  The  above  references,  therefore, 
do  not  consider  noise  and  are  directed  primarily  to  the 
calculation  of  signal  amplitude  at  the  output  of  the  limiter. 

The  purpose  of  this  paper  is  to  extend  the  Fourier- 
expansion  method  to  include  random  noise  in  the  analysis, 
and  to  derive  a  general  analytical  expression  for  the  output 
of  a  limiting  device.  The  approach  is  similar  to  that  used 
by  Reed  [9]  for  calculating  the  amplitudes  of  two  signals 
in  noise.  Closed-form  expressions  for  the  desired  signal 
and  cross-product  amplitudes  are  presented  for  the  case 
of  three  modulated  sinusoids  and  noise.  The  analysis  is  then 
extended  to  n  modulated  sinusoids  plus  noise,  and  approx¬ 
imate  expressions  for  the  signal,  cross-product,  and  noise 
terms  are  derived. 

Calculation  of  the  Limiter  Output 

The  specific  model  for  the  bandpass  limiter  to  be  con¬ 
sidered  is  shown  in  Fig.  1.  The  input  to  the  limiter 

x(t)  =  s(t )  +  n{t)  (1) 

consists  of  the  signal  s(t)  and  a  band  of  zero-mean  stationary 
Gaussian  noise  n(t).  It  is  assumed  that  the  bandpass  filter 
preceding  the  limiter  is  wide  enough  to  pass  the  signal  with 
negligible  distortion  and  limits  the  input  noise  to  a  narrow 
bandwidth  that  is  small  compared  to  the  center  frequency 
of  the  filter.  The  limiter  is  followed  by  another  bandpass 
filter  that  confines  the  output  spectrum  essentially  only 
to  the  fundamental  band  of  the  signal. 

It  is  assumed  that  the  limiter  has  a  hard-limiting  character- 
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Fig.  1.  Model  for  investigation  of  ideal  symmetric  limiting. 


istic  that  limits  its  output  to  either  ±1.  Thus,  if  the  limiter 
input  is  x(0,  the  output  y(t)  may  be  expressed  in  analytical 
form  [7]  as 

y(t)  =  -  f  sin  [tfx(0]  —  .  (2) 

'  n  J0  v 

When  the  expression  (1)  is  inserted  into  this  integral,  the 
sine  of  a  sum  may  be  expanded  into  the  sum  of  two  prod- 
ucts  of  sine  and  cosine.  Thus, 

9  r°°  dv 

y(t )  =  -  sin  [rs(0]  cos  IX(0]  — 

n  J0  v 

9  Z*00  dv 

+  -  cos  [as(0]  sin  [un(f)]  —  .  (3) 
nJ  o  ® 

The  narrowband  Gaussian  noise  at  the  limiter  input  may 
be  expressed  as 

n(t)  =  r  cos  ( co0t  +  <t>),  (4) 

where  the  envelope  r  and  the  phase  4>  are  slowly  time- 
varying  random  variables  having  Rayleigh  and  uniform 
distributions,  respectively. 

Substituting  the  above  expression  for  the  noise  in  (3), 
and  using  the  general  relationships 

sin  (z  cos  p)  =  2  £  (-l)mJ2m+1(z)cos(2m  +  l)p 

m  -  0 

cos  (z  cos  p)  =  JoO)  +  2  £  (-l)mJ2m(z)  cos  2 mp,  (5) 

m  -  1 

the  following  expression  for  the  limiter  output  is  obtained : 
9  r°°  dv 

y(t)  =  -  sin  0(0]  •  J0(vr )  — 
n  J0  v 

+  -  £  (-O'"  f  sin  0(0] 

7t  m  =  1  J  o 

•  J2Jvr)  —  cos  2m(o)0t  +  4>) 

v 

+  -  £  (-Om  f  cos  0(0] 

71  m  =  0  J  o 

•  J2m+i(yr)  —  cos  (2m  +  l)(oV  +  </>)•  (6) 

V 

All  the  terms  in  the  above  equation  are  random  functions 
due  to  the  presence  of  the  noise  envelope  r  and  phase  (j). 
However,  only  the  first  integral  will  yield  an  average  output, 
while  the  two  other  integrals  will  not  contribute  since  all 
the  terms  contain  the  random  phase  of  the  noise. 

The  average  limiter  output  at  any  arbitrary  time  t  is 
obtained  by  averaging  over  all  possible  noise  amplitudes, 


as  follows: 

Z(t)  =  £[y(0]  =  ;  f 
71  J0 


dv 


sin  [us(0]  '  £[70(rr)]  —  .  (7) 


Since  the  noise  envelope  has  a  Rayleigh  distribution 


p(r)  =  -=  exp 


2<r2J  ’ 


(8) 


where  <r2  is  the  total  noise  power  at  the  limiter  input,  the 
expected  value  of  Jo(vr)  is  given  by  [8] 


E[-Jo(”r)]  =  J  T0(t>r)p(r)  dr  =  exp  - 


(9) 


Substitution  of  (9)  in  (7)  yields 


This  is  identical  to  the  expression  obtained  by  Reed  [9] 
for  calculating  the  amplitudes  of  two  signals  in  noise.  The 
above  integral  may  also  be  written  in  terms  of  the  error 
function  [10]  as 


z(0  =  erf 


(U) 


The  presence  of  noise  at  the  input  of  a  hard  limiter  thus 
tends  to  make  the  limiting  soft,  due  to  the  gradual  saturation 
characteristics  of  the  error  function.  It  may  be  seen,  as 
also  observed  by  Jones  [2],  that  had  a  smooth  limiter 
having  an  error  function  amplitude  characteristic  been 
used,  the  only  change  in  (10)  and  (11)  would  have  been 
to  replace  <r2  by  cr2  +  y2.  The  quantity  y  determines  the 
slope  of  the  limiter  characteristic,  i.e.,  how  fast  saturation 
is  approached  in  the  case  of  soft  limiting.  The  effect  of  soft 
limiting  on  the  average  output  is  therefore  mathematically 
identical  to  adding  noise  to  a  hard  limiter,  provided  that 
an  error-function  representation  can  be  used  for  the  limiter 
characteristic. 

The  components  that  collectively  make  up  the  limiter 
output  noise  spectrum  are  the  second  and  third  integrals  of 
(6): 

t?(0  =  -  S  (-i)m  f  sin  1X0] 

7T  m=  1  Jq 


*  Jim(vr)  —  cos  2m((o0t  +  (j) ) 

V 


+  -  X  (-l)m  f  cos  [us(0] 

71  m  —  0  J  o 

•  Jim  +  i(ur)  — cos(2m  +  lXcOo*  +  <t>)-  (12) 

V 

The  expressions  (10)  and  (12)  are  the  two  fundamental 
equations  for  the  analysis  of  the  limiter  in  the  time  domain. 
In  Davenport’s  notation  [1],  (10)  represents  the  output 
signal  and  intermodulation  components  resulting  from  the 
interaction  of  the  signal  with  itself  (s  x  s  terms).  Similarly, 
(12)  represents  output  noise  components  that  are  due  to 
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the  interaction  of  the  input  noise  with  itself  (n  x  n  terms) 
and  with  the  input  signal  (s  x  n  terms).  A  separation  of  the 
noise  output  into  n  x  n  terms  and  s  x  n  terms  is  not 
possible,  without  the  knowledge  of  the  signal  waveform 
s(t). 

If  the  characteristic-function  method  had  been  used  to 
calculate  the  limiter  output,  the  phases  of  all  the  components 
at  the  limiter  output  would  have  been  lost,  and  a  separation 
of  the  output  components  as  in  (10)  and  (12)  would  have 
been  possible  only  in  terms  of  the  autocorrelation  function 
of  the  individual  components. 


where  hmt  represents  the  integral 


4  p  J,(va)Jn(vb)J((vc)  exp 
x  Jo  v 


t>y 
2  . 


dv.  (16) 


The  output  of  the  limiter  can  be  arranged  in  terms  of 
signal  and  cross-product  components  by  representing  each 
of  the  products  of  the  three  cosines  in  (15)  as  a  sum  of  four 
cosines  using  the  trigonometric  identity: 


cos  a  •  cos  ■  cos  y  =  ^[cos  (a  +  /J  +  y) 


+  cos  (a  +  p  -  y) 


Limiting  of  Three  Signals  and  Noise 

The  effect  of  ideal  bandpass  limiting  upon  one  and  two 
sinusoids  lying  in  narrowband  Gaussian  noise  has  been 
treated  by  several  authors.  This  section  presents  analytic 
results  for  the  case  of  three  angle-modulated  sinusoidal 
input  signals  plus  noise.  Let  the  signals  at  the  input  to  the 
limiter  be 


s(t)  =  a  cos  \co2t  +  <^(0]  +  b  cos  [ co2t  +  4>2(t) ] 

+  c  cos  [co3r  +  03(r)] 

=  a  cos  r  +  b  cos  s  +  c  cos  t,  (13) 

where  r  =  co,f  +  and  similarly  for  j  and  t. 

When  (13)  is  inserted  into  (10),  the  sine  of  a  sum  may  be 
transformed  by  simple  trigonometry  into  the  sum  of  four 
products  of  sines.  Thus 


v2a2  dv 
2  .  v 


[sin  ( va  cos  r)  cos  ( vb  cos  s )  cos  ( vc  cos  t) 

+  sin  (vb  cos  s)  cos  (vc  cos  t)  cos  (va  cos  r ) 

+  sin  (vc  cos  t)  cos  (va  cos  r )  cos  (vb  cos  s) 

—  sin  (va  cos  r )  sin  (vb  cos  s)  sin  (vc  cos  /)].  (14) 

Using  (5),  the  sines  and  cosines  may  be  expanded  as  Fourier 
series  in  r,  j,  t,  whose  coefficients  are  Bessel  functions.  The 
output  of  the  limiter  prior  to  filtering  may,  therefore,  be 
written  as 


zw  =  4i  z  i  (-i rm+<£m-£, 

»- 0  m  —  0  q- 0 


^2n+i,2mt2q  cos  (2w  +  l)r  cos  (2 ms)  cos  (2 qt) 


00  00  00 


+  4  I  I  I  (-l)n+m+«em  • 


'm  °<5j 


n  =  0  m  =  0  q  =  0 

h2q,2n+i,2m  cos  (2 qr)  cos  (In  +  l)s  cos  (2 mt) 

+  4  t  Z  f  (— l)"+m+,6m  •  e 

b=0  iti  —  0  9  =  0  ’ 

h2m,2q,2n+i  cos  (2 mr)  cos  (2 qs)  cos  (2n  +  1  )t 

4n?0  m?0  |0(-l)"  +  m  +  %B+1>2n+1.2,+  I 
cos  (2  n  +  l)rcos(2m  +  l)scos(2q  +  l)t,  (15) 


-I-  cos  (a  —  /?  +  y) 

+  cos  (a  -  p  -  y)].  (17) 

The  amplitude  of  any  desired  signal  or  cross-product 
component  may  be  obtained  from  (16).  The  signal  and  the 
strongest  cross-product  components  are  obtained  by 
setting  n  =  m  =  q  =  Oin  (15): 

z(0  *  h100  cos  r  +  h010  cos  s  +  h00i  cos  t 

—  Aul[cos  (r  +  s  +  t)  +  cos  (r  +  s  —  t) 

•  +  cos  (r  -  s  +  t)  +  cos  (r  -  s  -  *)].  (18) 

To  evaluate  the  coefficient  h^,  it  is  convenient  to  begin 
with  the  power-series  expansion  of  the  product  of  two  Bessel 
functions  [11]: 


Jq(vb)  ■  J((VC)  =  -  mmvcf  V  tM' 

'/•'  p=o  j>!(p  +  0! 

■  iFi  -P-Z,n  +  ,  (19) 

where  2Fi(  )  is  the  Gaussian  hypergeometric  function. 
Substitution  of  the  above  in  (16)  yields 


(zY  y  (-i y 

\2!  p=o  pi  (p  +  ^)! 


The  solution  of  the  above  integral,  which  is  attributed  to 
Weber  and  Sonine  [8],  is  given  in  terms  of  confluent  hyper¬ 
geometric  function.  The  coefficient  hmi  may  be  written 
explicitly  as 


hpqf  — 


sfesfehTfeb )'br.)' 


•  z 


(-1Y 


sr(,+ 


it  +  ij  +  ? 


) 


p=o  pi  (p  +  {)! 

■(±r ’M-*-'- 


ft  +  n  +  l  odd.  (21) 
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This  form  is  useful  for  numerical  computations  because 
the  Gaussian  hypergeometric  function  terminates  after  the 
(p  +  l)th  term  due  to  the  parameter  -p.  The  series 
representation  for  hmf  is  convergent,  unless  one  or  more 
of  the  input  signals  rise  substantially  above  the  input  noise 
level,  or  the  noise  is  weak  (<r2  «  0).1  In  this  case,  the  con¬ 
vergence  of  (21)  is  dependent  on  the  relative  amplitude 
of  the  three  signals.  In  the  general  case  of  three  signals  of 
arbitrary  amplitudes  in  weak  noise,  a  different  representa¬ 
tion  for  the  h  coefficients  is  required.  This  will  be  treated 
later  in  this  section. 

In  the  special  case  of  two  signals  having  equal  amplitude, 
(b  =  c),  the  expressions  for  may  be  simplified  by  using 
the  identity 


2Fi(a,P,y,l)  = 


r(y)T(y  -  «  -  ft) 
T(y  -  a)T(y  -  ft)’ 


y  #  0,  —  1,  —2,-  •  •, 


Re  (y  —  a  —  ft)  >  0  (22) 


.  2  /  a  Y(  b 

"*  7t(^!)  W2  J  W2  of  2 *1 

r  T(  — s) •  T(>?  +  |  +  2s  +  1) •  H*  +  (fi  +  a  +  01} 

’  JL  r(f/  +  s  +  i)  •  r(c  +  s  +  i)  •  r(»j  +  z  +  s  +  i) 


(25) 


For  large  values  of  signal-to-noise  ratio  a2/2cr2,  the 
asymptotic  expansion  of  (23)  for  the  confluent  hyper- 
geometric  function  can  be  used  for  evaluating  the  integral. 
Unfortunately,  the  integral  converges  only  for  the  first 
two  terms  of  the  series  expansion.  Thus,  the  expression  for 
h  t  will  be  valid  only  for  large  values  of  a2/ 2c1.  Using  the 
method  of  residues  for  solving  the  integral,  the  following 
expressions  for  the  amplitudes  of  the  desired  signal  com¬ 
ponents  are  obtained: 


for  the  Gaussian  hypergeometric  function.  If  the  input  noise 
can  be  neglected,  the  expressions  may  be  further  simplified 
by  using  the  asymptotic  expansion  for  the  confluent  hyper¬ 
geometric  function: 


,F i(«,  ft,  - Z ) 


m _ ‘ 

-  a)Za  . 


T(P  -  a)Z' 


i  ,  g(g  ~  ft  +  + 

Z 


Z  -y  00.  (23) 


_2_  s  rr(n  +  i)]3  M2n+1 

100  Ti7/2  „= o  (n  !)2(n  +  1) !  \2b) 

■  f-j—  +  3(i K«  +  1)  -  <A(n  +  i)}  +  21og  — 
In  +  1  a 

_2_/  ff  \2  f  fr (n  +  l)]3  /a\2B+1 
+  7t7/2V2b'  «=o  n !  [(n  +  1)  !]2  \2bl 


If  only  the  first  term  in  the  series  expansion  is  taken,  the 
expressions  for  the  amplitude  of  the  desired  signals  are  then 
identical  with  the  results  obtained  by  Sollfrey  [7].  Further¬ 
more,  we  also  obtain  expressions  for  the  amplitude  of  the 
cross-product  terms,  which  are  not  given  by  Sollfrey.  The 
expression  for  hmi  is  convergent,  however,  only  if  a  is 
greater  than  2b.  For  the  case  where  2b  is  greater  than  a  and  "oio 
noise  is  weak,  hmi  may  be  evaluated  by  the  method  used  by 
Sollfrey  for  the  noiseless  case.  It  involves  replacing  the 
product  of  the  two  Bessel  functions  of  the  same  argument 
in  (16)  by  a  contour  integral  representation 


Jn(pb)  ■  Jt(vb) 

-if 


r(-s)-rfa  +  j  +  2sjM) 


2ni  jl  T(f/  +  s  +  1)  •  r(S  +  s  +  1) '  J~(tj  +  £  +  s  +  1) 


(!) 


ij  +  $+2j 


ds 


(24) 


(26) 


•  [— - — - —  +  3{i j/(n  +  1)  -  +  i)} 

In  +  1  2n  +  1 

l  _  _2_  s  rr(n  +  Ul2  •  r(n  +  I)  (±\ 

n2  nV2  n=o  »![(«  +  I)!]2 


.( n  +  l)(2n  +  1) 

-  ift(«  +  i)}  +  2  log  —  I 

[-  2 


+  3{\p(n  +  1) 


2n  +  1 


[r(n  +  i)]2-r(n  +  I)M5 

(n  !)3  \2bJ 


+  3{t l/(n  +  1)  —  +  i)} 


and  then  solving  the  double  integral  by  changing  the  order 
of  integration,  i.e.,  by  evaluating  the  real  integral  over  v 
first.  The  real  integral  in  this  case  is  that  of  Weber  and 
Sonine  [8],  analogous  to  that  of  (20),  and  its  solution 
is  obtained  in  terms  of  the  confluent  hypergeometric 
function.  The  result  is 

1  The  series  in  (21)  is  convergent  for  all  parameter  values  except 

a2  =  0. 


+  2  log  —1  ,  ->1-  (2’). 

a  J  a 

Here,  )  denotes  the  logarithmic  derivative  of  the  gamma 
function.  Similar  expressions  can  also  be  derived  for  the 
desired  cross-product  components;  the  calculation  of  the 
residues,  however,  is  quite  tedious.  In  the  absence  of  input 
noise  (a2  =  0),  the  expressions  are  identical  with  the  results 
of  Sollfrey  [7]. 
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To  evaluate  the  coefficients  h ^  for  arbitrary  signal 
amplitudes  and  weak  noise,  it  is  convenient  to  use  Bate¬ 
man’s  expansion  [12]  for  the  product  of  two  Bessel 
functions : 


=  2bn  *  c*[2(b2  +  c2)]“  1/2(’+*+1) 


where 


•  £  f(ri,Z,n,b,c ) 

«”  0 

.  n+i(v^2(b2  +  c2)) 

i? 


tions.  For  example,  the  expression  for  hl00 , 
2  ' 


^100  “ 


2a2 


71  V  (h2  +  c2)  n^o 
•  2Fj  n  +  1,  1, 


I  (-l)n2Fi(-«,  «  +  1,  l,i) 
) 


b2  +  c2 

lF>  (n  +  i'-"-h2-WTs) 

may  be  simplified  by  using  the  identity 


(32) 


f{U,n,b,c)  =  Lt.1^  +A±2n±  1)  •  r(q  +  £  +  n  +  1)-IT({  +  »  +  1) 

ni-rfa  +  k  +  i)-r(£  +  i)-r({  +  i) 


•  2^i(-n.  >/  +  £  +  n+i,  £  +  1, 

Both  the  hypergeometric  functions  terminate  after  the 
(n  +  l)th  term  due  to  the  parameter  -n. 

Substitution  of  (28)  in  (16)  yields 


i)  •  iF i  r\  +  £  +  n  +  1,  £  +  1,  j)  .  (29) 

»  +  i,  i, i)  “  '/”  [r  (|  -  !•)  ■ r  (i  +  *)]"'• 


*«  p (I,*  +  6,0 

7T  n  =  0 

.  f°°  Jn(Va)  *  ^ij  +  g+2n  +  j(V\/2(b2  +  C2)) 

Jo  r2 

[i;2cr2l 

"  (30) 

The  integral  can  be  easily  evaluated  by  replacing  the 
product  of  the  Bessel  functions  by  the  series  representation 
of  (19).  In  the  absence  of  noise,  the  integral  represents  the 
discontinuous  integral  of  Weber  and  Schafheitlin  [13]. 
Its  solution  depends  on  the  relative  magnitude  of  the 
arguments  of  the  two  Bessel  functions.  In  the  following, 
only  the  solution  for  a  <  \/2 (b2  +  c2)  will  be  considered; 
the  case  a  >  y]2(b2  +  c2)  can  be  treated  similarly.  Using 
the  appropriate  solution  of  the  integral  [13],  the  expression 
for  h ^  may  be  written  as 

Knt  =  ;  •  »  ■  cW2  +  c2)]' 1/2("+,’+«  £  f(U,n,b,c) 

IE n  +  JQ 1+JL+  01 

n\  •  r[n  +  2  +  i(t}  +  £  -  p)] 

•  iFi  +  i(fi  +  rj  +  £),  -n  -  1  +  i(/i  -  n  ~  £)> 

M  +  1  w+  c2))  ’  V2(fc2  +  c2)  >  a-  (31) 

This  expression  is  convergent  and  yields  the  amplitudes 
of  the  signals  and  the  cross  products  in  the  absence  of  noise. 
In  most  cases,  it  would  be  possible  to  simplify  the  expression 
for  any  h  coefficient  by  expressing  the  hypergeometric 
function  containing  the  factor  \  in  terms  of  gamma  func- 


(33) 

It  is  seen  that  the  series  is  zero  for  all  odd  integer  values  of 
72.  In  the  special  case  of  two  signals  having  equal  amplitudes, 
the  second  hypergeometric  function  can  also  be  replaced 
by  (33),  yielding  the  simple  expression 


*100 


1  (“)  y  1T(*  +  i)T 

n2  \b)  n—0  (n!)2 


—  2tj  —  2 


)■ 


-  >  1.  (34) 


Similarly,  for  h0l0  the  following  expression  is  obtained 


^oio  —  h 


ooi 


1  f  (-i)n-r(72  +  i ) 

7I2  n=o  (n  +  1)  *  T(n  +  1) 


a 


by  using  the  expression  in  the  brackets  to  replace 
W n  •  (n  +  1)  *  2F\ (  —  «,  72  +  2,  1,  £).  Numerical  evaluation 
of  (34)  and  (35)  yield  results  identical  to  those  obtained 
from  (26)  and  (27)  in  the  absence  of  noise;  however,  it  has 
not  been  possible  to  convert  one  form  into  the  other. 
Similar  expressions  for  the  cross  products  can  be  readily 
derived.  This  approach  has  the  advantage  over  Sollfrey’s 
method  in  that  it  directly  provides  the  expressions  for  the 
signals  and  all  the  cross  products.  Thus,  there  is  no  need  to 
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calculate  the  residues  for  each  h  coefficient,  as  required 
in  the  evaluation  of  (25). 

Numerical  Results 

The  results  for  limiting  of  three  signals  in  noise  are  shown 
in  Figs.  2-6.  In  all  cases,  the  amplitudes  of  two  signals  were 
assumed  to  be  equal  (b  —  c ),  in  order  to  reduce  the  number 
of  variables.  The  following  expressions  were  used  for 
numerical  computations: 


/  a' 

1  y  (" 

iy-(2p)\-T(p  +  i) 

l  b 

:  W2 

(p'.f 

W2  i 

J 

(p 

+  2, 

2c1) 

(36) 

■(  b  ' 

1  y  (= 

i)p-(2  P  +  i)!-r(p  +  i) 

:  W2  O' 

1  p = 0 

Q>!)2-[(p  +  l)!]2 

(V2.) 

2  P 

1^1 

(P  +  h  ''  -  2?)  ■ 

(37) 

These  are  obtained  from  (21)  by  replacing  the  Gaussian 
hyper  geo  metric  function  with  the  identity  of  (22).  When 
the  noise  is  weak,  computation  of  these  expressions  breaks 
down  and  it  is  necessary  to  employ  the  asymptotic  expan¬ 
sion  of  (23)  for  the  confluent  hypergeometric  function. 
The  resulting  expressions,  however,  are  convergent  only 
if  a  is  greater,  than  2b.  For  the  case  a  less  than  2b  (26)  and 
(27)  are  applicable. 

The  curves  in  Figs.  2  and  3  show  the  amplitudes  hl00 
and  h010  plotted  against  the  limiter  input  power  ratio  of 
double  to  single  component  with  the  limiter  input  signal- 
to-noise  ratio  of  the  single  component  as  parameter.  For 
bja  small,  that  is,  strong  single  and  weak  double  component, 
the  results  approach  the  one-signal-in-noise  case  discussed 
by  Davenport  [1].  In  the  absence  of  noise,  the  single  com¬ 
ponent  tends  to  4/7i  and  gets  essentially  all  of  the  limiter 
output  power.  As  the  signal-to-noise  ratio  decreases,  the 
noise  uses  up  some  of  the  output  power,  so  that  the  pro¬ 
portion  used  by  the  strong  single  component  decreases. 
As  the  ratio  bja  increases,  the  proportion  of  the  output 
power  used  by  the  double  component  also  increases,  while 
that  of  the  single  component  decreases.  When  the  double 
component  is  strong,  that  is  bja  is  large,  the  results  approach 
the  two-signal  case  discussed  by  Jones  [2].  The  double 
component  tends  to  8/7i2,  when  there  is  no  noise.  Each 
of  the  two  strong  signals  receives  —3.9  dB  of  the  output 
power,  so  that  0.9  dB  remains  for  the  weak  single  com¬ 
ponent  and  cross  products.  As  the  signa!-to-noise  ratio 
decreases,  the  share  of  the  output  power  consumed  by 
noise  increases,  and  consequently  the  power  in  the  double 
component  decreases.  The  plot  in  Fig.  4  shows  the  result 
of  limiting  three  signals  of  equal  amplitude  in  noise. 

Signal  suppression  is  another  phenomena  of  interest. 
It  is  defined  as  the  ratio  of  the  weak-to-strong  component 
ratio  at  the  limiter  output  to  the  corresponding  quantity 


Fig.  2.  Amplitude  of  the  limiter  output  single  component — three 

inputs. 


Fig.  3.  Amplitude  of  the  limiter  output  double  component — three 

inputs. 


at  the  input.  In  decibels  the  suppression  ratio  is  given  by 


L  =  20  log  hi0° 

-  20  log  l , 

-  >  i 

^010 

b 

b 

=  20  log  ^ 

-  20  log  -  , 

-  <  1. 

(38) 

^100 

a 

b 

Positive  decibel  values  of  L  correspond  to  reduction  of  the 
weak  component.  The  first  case  corresponds  to  a  strong 
single  component  and  weak  double  components,  the 
second  to  weak  single  and  strong  double.  In  Figs.  5  and  6, 
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Fig.  4.  Amplitude  of  the  limiter  output  signal  components,  three 
inputs  equal  in  amplitude. 


the  suppression  ratio  L  is  plotted  against  the  input  ratio 
bja  in  absolute  decibels. 

For  the  one-strong-two-weak  case  shown  in  Fig.  5,  the 
maximum  suppression  of  the  weak  double  component 
occurs  at  large  input  signal-to-noise  ratios,  and  it  is  seen 
to  be  6  dB  when  there  is  no  noise.  On  the  other  hand, 
when  the  noise  is  strong,  the  limiter  acts  essentially  as  a 
linear  device,  in  the  sense  that  almost  no  suppression  of  the 
weak  double  component  takes  place.  This  behavior  is  in 
complete  agreement  with  the  suppression  of  a  weak  signal 


by  a  strong  signal,  as  influenced  by  the  noise  level,  for  the 
two-signal  case  discussed  by  Jones  [2].  In  the  case  of  strong 
double  and  weak  single  components,  shown  in  Fig.  6, 
hard-limiting  enhances  the  weaker  component  with  respect 
to  the  stronger  double  components  at  large  values  of  bja. 
Thus,  under  these  circumstances,  the  limiter  displays 
“negative  suppression.”  Maximum  enhancement  of  the 
weak  signal  occurs  in  the  absence  of  noise. 

This  interesting  phenomenon  was  originally  discovered 
by  Jones  [2]  for  the  case  where  the  weak  signal  was  noise 
rather  than  another  sinusoid.  In  the  case  of  three  signals, 
the  effect  arises  from  beats  between  the  two  equal  strong 
components,  which  permit  the  weak  component  to  slip 
through.  It  also  appears  in  the  results  of  Shaft  [3]. 

The  asymptotic  curves  in  Figs.  2-6  describing  the  noise¬ 
less  case  were  first  obtained  by  Sollfrey  [7]. 


Limiting  of  n  Signals  and  Noise 

The  method  used  to  analyze  the  case  of  three  input  signals 
may  be  extended  to  the  general  case  of  n  angle-modulated 
sinusoids.  The  signal  and  intermodulation  products  at  the 
limiter  output  are  obtained  by  inserting  the  expression  for 
the  limiter  input  signals 


n 


s(t)  =  £  Oi  COS  [ COit  +  0,(0] 

i  =  1 


(39) 


into  (10)  and  using  the  general  relationship 
sin  ft  cos 


=  I 


•  E 

Pn=  -  00 


n  sin  (i  Pi  (<*.-  +  . 


(40) 

where  the  symbol  J“I  indicates  that  all  n  of  the  JPi  coeffi¬ 
cients  are  multiplied  together.  Thus,  the  limiter  output, 
prior  to  bandpass  filtering,  is  given  by 


*(')-*  £  -  £ 

Pi-  -00  pn=-00 


•  sin 


(Oit  +  <t>i{t) 


where 


h 


PlP2*  •  * Pn 


JPi(vai)  exp 


(42) 

2  J  v 


represents  the  amplitude  coefficient  of  the  output  signal 
and  intermodulation  components. 

It  is  apparent  that  the  determination  of  all  those  terms  in 
(41)  that  will  fall  in  the  passband  of  the  filter  is  quite  a 
formidable  task  when  n  is  large.  Fortunately,  Shaft  [3] 
has  shown  that  the  dominant  output  components  are  the 
signal  and  cross-product  terms  for  which  is  small 

(<  2).  This  approximation  considerably  simplifies  (41), 
since  instead  of  considering  all  the  pt  from  ~oo  to  oo,  it 
is  only  necessary  to  consider  them  in  the  range  from  —2 
to  2.  The  significant  filter  output  components  may  be 
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determined  by  first  considering  the  different  frequency 
combinations  in  (41),  and  then  selecting  the  desired  ones, 
i.e.,  the  combinations  that  will  yield  frequencies  in  the 
passband. 

A  closed-form  solution  for  the  amplitudes  of  the  output 
components  is  extremely  difficult  when  n  is  large.  For  more 
than  three  signals,  the  amplitudes  may  be  determined 
either  by  numerically  evaluating  the  integral  of  (42),  or  by 
using  approximations  to  solve  the  integral.  Shaft  [3]  used 
the  first  approach  and  has  reported  extensive  computer 
results  for  as  many  as  a  hundred  input  signals.  Approx¬ 
imate  expressions  for  the  amplitudes  of  the  signals  and 
cross  products  may  be  obtained  by  using  approximations 
for  the  Bessel  functions.  As  an  illustration,  consider  the 
signal  at  frequency  f1.  The  amplitude  and  phase  of  the 
signal  is  obtained  from  (41)  and  (42)  by  setting  py  =  ±  1 
and  all  the  other  pt  equal  to  zero.  The  result  is 

zx(t)  =  hio- ■  -o  cos  (coj?  +  <t>i(t)),  (43) 


where 

*10 


....  -ir 

n  J  o 


Jlivai)  J]  Jo(ufli)  exp - — 


[-f: 


dv 

v 


(44) 


An  approximate  solution  to  (44)  is  obtained  by  using  an 
exponential  approximation  for  the  product  of  Bessel 
functions,  as  suggested  by  Gyi  [4] : 


fl  J0(va,)  «  exp  [-i^i2], 

i  -  2 


where 


The  result  is  then 


(45) 


(46) 


where  p  =  a^i\l{G2  +  c^2)].  Bounds  on  the  error  for 
this  approximate  solution  of  Gyi  have  been  recently  ob¬ 
tained  by  Campbell  [15].  The  bound  is  of  the  order  n"1 
when  the  amplitudes  of  all  the  signals  are  approximately 
equal. 

The  strongest  cross-product  terms  are  those  that  are 
produced  due  to  the  mixing  of  any  two  or  three  input  signals. 
Approximate  expression  for  the  amplitude  of  these  cross 
products  can  be  obtained  from  the  results  for  two  signals 
and  three  signals,  respectively.  As  an  example,  consider 
the  strongest  cross  products  due  to  the  mixing  of  the 
frequencies  fl9  f2i  and  /3.  They  may  be  obtained  from 
(41)  and  (42)  by  considering  all  the  combinations  of  these 
three  frequencies.  The  result  is 

3W(0  =  -Amo  -  o[c°s  (r  +  s  +  t)  +  cos  (r  +  s  -  t) 


+  cos  (r  -  s  +  t)  +  cos  (r  -  s  -  f)]>  (47) 


where  r  =  u>xt  +  and  similarly  for  s  and  t.  The 

amplitude  coefficient  is  given  by 

j  /*  oc  n 

*1110- - -o  =  -  Jl(va1)Jl(va2)J1{vai)  fl  7o(^,) 
n  J  >=4 


v2a2' 


dv 

v 


.  (48) 


Using  the  exponential  approximation  to  the  product  of  the 
Bessel  functions,  the  solution  of  the  above  integral  may  be 
obtained  from  (21)  merely  by  replacing  a2  by  a2,  where 


The  components  that  collectively  make  up  the  limiter 
output  noise  may  be  obtained  by  substituting  (39)  for 
s(t)  in  (12),  and  expanding  sin  [ttf(0]  and  cos  [tw(0]  with 
the  aid  of  (40).  The  expansion  of  cos  [t>s(f)]  contain 
cosine  terms  instead  of  sine  in  (40).  The  resulting  expression 
then  represents  the  combined  output  noise  produced  due 
to  the  interaction  of  the  input  noise  with  itself  («  x  n 
terms)  and  that  produced  as  a  result  of  the  interaction  of 
the  input  signal  with  noise  (j  x  n  terms).  In  this  expression, 
the  term  corresponding  to  m  and  all  pl  equal  to  zero 
represents  the  direct  feedthrough  noise  (n  x  n  terms),  and 
all  the  remaining  terms  constitute  the  noise  produced  by 
the  mixing  of  the  signal  with  the  limiter  input  noise  (s  x  n 
terms).  The  expression  for  the  direct  feedthrough  noise 
(n  x  n  terms)  is  given  by 


a  r°°  ~  dv 

rinx„(t)  =  -  Ji(»r)  n  Jo(vad  —  (cos  co0t  +  <f>).  (49) 
n  J0  i=i  v 

Davenport  [1]  has  shown,  for  the  case  of  one  signal  and 
noise,  that  only  the  n  x  n  terms  contribute  significantly 
to  the  output  noise  at  low  limiter  input  signal-to-noise 
ratios.  When  the  number  of  signals  is  sufficiently  large, 
so  that  the  amplitude  distribution  of  their  sum  is  approx¬ 
imately  Gaussian,  the  signal-to-noise  ratio  for  any  one 
signal  will  be  small  at  the  limiter  input.  Consequently, 
from  Davenport’s  results,  it  is  reasonable  to  expect  that 
only  the  direct  feedthrough  noise  will  contribute  signifi¬ 
cantly  to  the  limiter  output  noise. 

Again,  using  the  exponential  approximation  to  the 
product  of  the  Bessel  functions  (49)  can  be  readily  integrated 
[8]  to  yield 

-  Jl  (f )  'F'  (*• 2-  -  0) cos("0' +  m' 

(50) 


where 


represents  the  total  power  of  all  the  input  signals. 

The  power  in  the  feedthrough  noise  component  is  given 


>340 


IEEE  TRANSACTIONS  ON  INFORMATION  THEORY,  MAY  1972 


by 

JV0  =  ~  E  [  f  J^vr)  exp  [-^2cr32]  — 

7T2  L  J0  1>J 

= 4  r  r exp  + v2n 

Tt  JO  Jo 

■  E[J i(ur)  •  J jO*)]  —  —  ,  (51) 

U  V 

where 

E[Ji(ur)  ■  J^vry]  =  —  f  J,(ur)  •  Jt(vr)  ■  r 
a  Jo 

'exp[~  £~]dr' (52) 

This  integral  is  well  known  as  Weber’s  second  exponential 
integral  [14],  Using  its  solution  [14],  (51)  may  be  written 
explicitly  as  [16] 


a  time-domain  approach.  A  major  advantage  of  the  time- 
domain  method  is  that  it  preserves  the  phases  of  the  signal 
and  intermodulation  products,  which  are  destroyed  in  the 
characteristic-function  method  generally  used  in  limiter 
studies.  Expressions  for  the  desired  signal  and  inter¬ 
modulation  product  amplitudes  have  been  obtained  for  the 
case  when  the  input  consists  of  three  angle-modulated 
sinusoids  and  noise.  The  analysis  is  extended  to  n  angle- 
modulated  sinusoids  plus  noise,  and  approximate  expres¬ 
sions  for  the  signal,  intermodulation  product,  and  noise 
terms  are  derived.  It  has  been  possible  to  obtain  a  time- 
domain  representation  for  the  limiter  output  noise,  which  is 
particularly  useful  for  the  evaluation  of  the  system  following 
the  bandpass  limiter. 
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=  1  f  [r(n  +  l)]2  /  <r2  \2n+1 

7T2  n  =  o  n  \  •  (n  -F  1)!  \(72  +  cr32/ 


(53) 


It  is  seen  that  (50)  represents  the  dominant  output  noise 
in  the  time  domain.  This  representation  is  particularly 
useful  for  the  evaluation  of  the  system  following  the  band¬ 
pass  limiter.  Equation  (50)  may  also  be  expressed  in  the 
form 


tlnxniO  =  *i(0  cos  co0t  -  yx(t)  sin  c o0t,  (54) 

where  xA(/)  and  yx(t)9  the  amplitudes  of  the  “in-phase” 
and  “quadrature”  components,  are  given  by 


The  functions  ^x(/)  and  >>i(0  are  in  general  non-Gaussian 
random  variables,  whose  distributions  may  be  obtained 
from  the  distribution  of  r  and  4 >• 

Conclusions 

The  effect  of  ideal  bandpass  limiting  on  an  arbitrary 
signal  lying  in  narrowband  stationary  zero-mean  Gaussian 
noise  has  been  analyzed,  and  general  analytic  expressions 
for  the  limiter  output  components  have  been  derived  using 
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